Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN WF$t style file v2.2) 



Detection of the radial velocity curve of the B5— AO 
supergiant companion star of Cir X— 1? 

P.G. Jonker 1 ' 2 ' 3 *, G. Nelemans 4 , C.G. Bassa 3 

- - - 1 SRON, Netherlands Institute for Space Research, Sorbonnelaan 2, 3584 CA, Utrecht, The Netherlands 

i 2 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, Massachusetts, U.S.A. 

• 3 Astronomical Institute, Utrecht University, P.O.Box 80000, 3508 TA, Utrecht, The Netherlands 

' 4 Department of Astrophysics, IMAPP, Radboud University Nijmegen, Toernooiveld 1, 6525 ED, Nijmegen, The Netherlands 

CN ■ 
-i— > 
O 

O 



5 February 2008 



ABSTRACT 

CNJ ■ In this Paper we report on phase resolved /-band optical spectroscopic and photomet- 

ric observations of Cir X— 1 obtained with the Very Large Telescope. The spectra are 
dominated by Paschen absorption lines at nearly all orbital phases except near phase 
zero (coinciding with the X-ray dip) when the absorption lines are filled-in by broad 
Paschen emission lines. The radial velocity curve of the absorption lines corresponds 
to an eccentric orbit (e = 0.45) whose period and time of periastron passage are con- 
' sistent with the period and phase predicted by the most recent X-ray dip ephemeris. 

f*^ , We found that the J-band magnitude decreases from 17.6 to ~ 16.8 near phase 0.9- 

1.0, this brightening coincides in phase with the X-ray dip. Even though it is likely 
that the absorption line spectrum is associated with the companion star of Cir X-l, 
we cannot exclude the possibility that the spectrum originates in the accretion disc. 
However, if the spectrum belongs to the companion star, it must be a supergiant of 
Jjj | spectral type B5-A0. If we assume that the compact object does not move through 

the companion star at periastron, the companion star mass is constrained to < 10 Mq 
for a 1.4 M© neutron star, whereas the inclination has to be > 13.7°. Alternatively, the 
measured absorption lines and their radial velocity curve can be associated with the 
accretion disc surrounding a 1.4 Mq neutron star and its motion around the centre 
of mass. An absorption line spectrum from an accretion disc is typically found when 
our line-of-sight passes through the accretion disc rim implying a high inclination. In 
this scenario the companion star mass is found to be ^0.4 Mq. However, from radio 
observations it was found that the angle between the line-of-sight and the jet axis is 
smaller than 5°. This would mean that the jet ploughs through the accretion disc in 
this scenario, making this solution less probable. 

Key words: stars: individual (Cir X-l) — accretion: accretion discs — stars: binaries 
— stars: neutron — X-rays: binaries 
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1 INTRODUCTION 

X-ray binaries are binary systems harbouring a neutron star 
or black hole compact object that accretes matter from ei- 
ther a low- or a high-mass companion star (LMXB and 
HMXB, respectively). LMXBs are typically old systems 
whereas the early type companion star of HMXBs precludes 
systems older than a few times 10 7 yr. Recently, it has 
been realised that many X-ray binaries might have started- 
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off as intermediate mass X-ray binaries (cf. C yg X-2; 
Tau ris fc Savoniielll999tlPodsiadlowski et al]l2002l) . One in- 
triguing system that so far has defie d classification is Cir X- 
1. I t has a 16.6 day orbit al period llKaluzienski et af]ll97& 
SCI. 1 IClarkson etaH 120041 for the latest X-ray ephemeris). 
Owing to the detection of type I X-ray burst s, the compact 
object in Cir X-l likely is a neutron star jTennant et alJ 
ll986bl:lTennant et aljll986al) . 

The X-ray and radio behaviour of the source is complex. 
The long-term X-r ay lightcurve of the source has been dis- 
cussed in detail bv lSaz Parkinson et alJ <l2003l) . One of the 
striking features of the X-ray lightcurve is the periodic ap- 
pearance of dips. In radio, an arc minute scale radio nebula 
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was found dStewart et al.lll993ft . Furthermore, a relativistic 
outflow has been detected on scales of arc seconds that is 
aligned with the a rc minute scale jet ijFender et, al.l fl998t 
iFender et aljl2004T) . The inclination of the jet with respect 
to th e line-of-sight has to be less than 5° llFender et alJ 
2004). The presence of relativistic outflows detectable as 
synchrotron emission in the radio band ar gues against a high 
magnetic field neutron star in Cir X-l fcf. lFender fc Hendrvl 
2000). Furthermore, the detection of type I X-ray bursts 
from Cir X-l suggests that the neutron star has a low mag- 
netic field, since the thermonuclear instability giving rise to 
type I X-ray bursts is suppressed by dipole magnetic fields 
£ 10 12 Gauss Jjoss fc Lll980llBildstenll995l) . In accordance 
with this, type I X-ray bursts are not found in HMXBs. The 
low magnetic field of the neutron star in Cir X-l suggests 
that it is old and hence the companion star is not an early 
type star. 

So far, clear spectroscopic evidence on the nature of 
the companion star is lacking. It has been suggested tha t 
Cir X-l has a supergiant comp anion jMurdin et alJll98Cl) . 
However, it can be argued (cf. IJohnston^tHdT^)9£ ) that 
the source is too faint in the optical bands for the distance 
of 8-10 kpc that has been deriv ed from the t ype I X-ray 
bursts \ Jonker fc Nelemansl2004T) . Nevertheless. IClark et all 
( 2003) found emission lines in the near-infrared spectrum 
that are consistent with a mid-B super giant. However, these 
authors concluded that these features must have arisen in 
the accretion disc and/or in outflows from the disc. 

In this Manuscript we present phase resolved Very Large 
Telescope (VLT) photometric and spectroscopic observa- 
tions obtained with the FOcal Reducer/low dispersion Spec- 
trograph 2 (FORS2). 



2 OBSERVATIONS, ANALYSIS AND RESULTS 

We have obtained VLT/FORS2 spectra of the peculiar X- 
ray binary Cir X-l using the 1028z holographic grism with 
a slit width of 1". The observations were obtained in service 
mode on 21 different nights in the period ranging March 15- 
May 15, 2005 (MJD 53446-53507). In order to sample the 
~16.6 day long binary orbital period of Cir X-l we obtained 
one spectrum per night with an exposure time of 1730 sec- 
onds. Exceptions were May 13 and April 4. On May 13 we 
obtained 9 spectra with an exposure time of 1675 seconds 
since that night Ci r X-l was close to pe riastron (according 
to the ephemeris of lClarkson et af]l2004f) . On April 4, three 
spectra were obtained since the seeing was higher than the 
specified conditions on 2 of the 3 occasions. To minimise the 
light coming from an unrelated nearby field star (star 2 in 
iMonetllQga) . observations were obtained under good seeing 
conditions (seeing between 0'.'4 and l'.'3 as measured from 
the point spread function full-width-at -half-maximum of 
the acquisition images). The dispersion was 0.86 A per pixel. 
With the 1" slit width the resolution is about 120 km s _1 
at 8800 A. 

The spectra have been reduced with iraf 1 . We used 
the overscan area of the Charge Coupled Device (CCD) for 

IRAF is distributed by the National Optical Astronomy Obser- 
vatories 



bias subtraction. The data we re flatfield corrected and op- 
timally extracted dHornell986l) . Wavelength calibration was 
done using lines from He, Ar & Ne lamp spectra that were 
obtained during daytime the day after the observations with 
the same instrument set-up, as is customary for VLT Ser- 
vice mode observations. The rms scatter of the wavelength 
calibration was in the range of 0.05-0.08 A. 

The extracted spectra were further reduced and anal- 
ysed using the software package molly. We corrected the 
wavelength calibration for potential shifts caused by flexure 
by cross correlating the spectra over the wavelength ranges 
9050-9160 A and 9300-9500 A with the first object spec- 
trum. That part of the spectrum is dominated by features 
from the night sky which should have the same wavelength 
in each spectrum. Note that this does not account for un- 
certainties in the wavelength calibration caused by centroid- 
ing or star tracking inaccuracies. We did this for both the 
spectra of Cir X-l and that of star 2 that was also in the 
slit. Next, the observation times were corrected to the He- 
liocentric Julian Date time frame (using UTC times) and we 
normalised and rebinned the spectra to a uniform velocity 
scale removing the Earth's velocity. Cross-correlation of the 
spectra of star 2 with the first spectrum of that star over the 
ranee 8000-9000 A shows that the remaining rms velocity 
differences are 2.5 km s -1 . 

In Fig. Q we show the normalised spectra as a func- 
tion of the orbital period. To fold the spectra we used our 
best-fit epoch of periastron (T, see below) and an orbital 
period of P or b = 16.54 days as found b y extra polating the 
X-ray dip ephemeris of IClarkson et all i2004h . We binned 
the data in 15 phase bins. However, none of the spectra fall 
in the phase range 0.60-0.85, hence only 11 phase bins are 
shown. The most striking feature is the large change in the 
profile of the Paschen lines. At phase zero, large, broad (pos- 
sibly double peaked) emission lines are present. Superposed 
is an absorption line spectrum that could originate in the 
companion star. This absorption line spectrum can be seen 
clearly at phases 0.2-0.6. Besides the Paschen absorption 
lines that are often seen in B and A stars, small absorption 
lines oft en seen in supergiants can be seen as well, e.g . near 
8777 A jMunari fc TomasellallT999l : IClark et alJl2005t com- 
pare the Cir X-l spect rum with that of the su pergiant sys- 
tem XTE J 1739-302 in lNeeueruela et all2006h . This feature 
is most likely due to He I. 

We investigated the behaviour of the equivalent width 
of the luminosity indicator Paschen-12 at 8750 A as a func- 
tion of the orbital phase. At phase 0.4-0.6 the equivalent 
width of the absorption line is largest. Assuming that the 
Paschen lines are formed in the companion star, the con- 
tamination of the accretion disc at phase 0.4-0.6 is thus 
minimal but not necessarily zero. Hence, the measured 
equiva lent width of 2 A at phase 0.4-0.6 is a lower limit . 
Using iDanks fc Dennefeldl Jl994h. lAllen fc Strorr] <199ElL 
ISilva fc Corned <1992ft and lMunari fc Tomasellal <1999|) we 
find that the spectra at phase 0.4-0.6 resemble most that of 
late B/early A stars. As can be seen in Fig. [5] the lines are 
very narrow and from this we conclude that the star must be 
a supergiant. However, we cannot exclude that the absorp- 
tion line spectrum is from an accretion disc. Nevertheless, in 
order to test what spectral type would describe the observed 
spectrum best under the assumption that the absorption 
lines are stellar, we obtained template stellar spectra from 
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Table 1. Template stars from lLe Borgne et alJ 



Spectral type & Henry Draper Identification 



BHb 
B5Ib 
B7Iab 
AOIb 



HD 091316 
HD 164353 
HD 268749 
HD 087737 



B2Ia 
B3Ia 
B9Ia 
A3Ia 
F5Ia 



HD 268623 
HD 271163 
HD 032034 
HD 033579 
HD 269697 



iLe Borgne et al] (120031) . We optimally subtracted the spec- 
tra of 9 stars (see Table Q with spectral type ranging from 
Bllb-F51a from the Cir X-l » 8250 - 9300 A spectrum at 
phase 0.4-0.6. The optimal subtraction is performed on nor- 
malised spectra, minimising the residuals using the following 
recipe: f^ lrX_1 — A + bf^ cmp , where A is the assumed accre- 
tion disc contribution and b is the fraction of light from the 
companion star. From the residuals of the optimal subtrac- 
tion it was clear that the la supergiants matched the narrow 
absorption lines better than those of the lb stars. Between 
the la supergiants we found that the spectra of those of spec- 
tral type mid-B provided the best fit. The multiplicative 
factor in the optimal subtraction was consistent with unity. 
However, systematic uncertainties prevent us from obtain- 
ing formal (reduced chi-squared, xt) fits. Since the template 
star spectra and the Cir X-l spectrum are not obtained with 
the same telescope-instrument combination the instrumen- 
tal broadening of the line profile is different. Furthermore, 
the broadening due to the rotational velocity of the stars is 
different for the different stars. 

Next, we investigated the velocity of the absorption line 
spectrum as a function of the binary orbital phase. Since we 
did not obtain early type supergiant template star spectra 
we cross correlate the spectra of Cir X-l with that of one 
of the spectra. In the cross correlation it is vital that the 
continuum level has been determined accurately. Further- 
more, one has to avoid regions of the spectrum affected by 
skylines. In order to avoid problems with the continuum nor- 
malisation and skylines we have cross-correlated the region 
between 8700-8900 A. In Fig.|3]we show the observed radial 
velocity curve. The solid line in Fig. [3]is the best-fit ellipti- 
cal orbit (see Tableland |3J. However, in the fit shown in 
Fig. |3]we have increased the errors by 6.5 km s _1 in order 
to obtain a xl of 1.0. As mentioned above the rms residual 
velocities measured by cross correlating the spectra of the 
nearby star with itself are 2.5 km s" 1 . We used this 2.5 km 
s" 1 as a measure of the amplitude of systematic effects. This 
systematic error dominates the statistical error of the cross- 
correlation of the Cir X-l spectra. However, even when we 
take this uncertainty into account, the formal ~ 6 for the 
25 degrees of freedom. Possible other systematic effects that 
are not included are: (i) the effects of differences in X-ray 
heating on the absorption line spectrum. E.g. it can be seen 
comparing the top and bottom left panel of Fig. |3]that the 
first three velocity measurements are somewhat below the 
fit, these observations were performed when the X-ray flux 
was still relatively high (ii) from Fig. Q it can be seen that 
near phase zero strong Paschen emission lines are present. 
These emission lines are redshifted with respect to their rest 
wavelengths and since the Paschen absorption lines at those 
orbital phases are red shifted by a different amount, the 
emission lines fill-in the absorption in an asymmetric way, 
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Figure 1. VLT/FORS2 spectra of Cir X-l ranging from 8300- 
8900 A for the different orbital phases as indicated on the right 
hand side in the plot. The Paschen lines are easily discernible. 
They change from absorption to strong emission lines (with an 
absorption core) near phase zero. An arbitrary offset has been ap- 
plied to the spectra for display purposes. The X-es mark regions 
where cosmic ray hits have been incompletely corrected for. The 
absorption line spectrum likely reveals the companion star. The 
line at 8619.5 A is due to a diffuse interstellar absorption band 
(DIB). The feature near 8650 A is also due to DIB although an 
He I line at 8648 A found in B stars may contribute to the broad 
nature of the line. 



possibly skewing the cross correlation velocities (iii) slit cen- 
troiding and tracking errors can have introduced errors in 
the wavelength calibration. 

Since we cross correlated the spectra with a spectrum of 
Cir X-l itself the systemic velocity that is derived from the 
elliptical orbit fit (3.7±1.3 km s _1 ) is not the true systemic 
velocity. In fact, as one would expect, within 3 a it is consis- 
tent with 0. The systemic velocity is ~ —26 ± 3 km s -1 . We 
derived this by comparing the best-fit central wavelength 
of the Gaussian absorption lines with rest wavelengths of 
8467.25, 8598.4, 8750.5 and 8862.8 A. The Gaussians were 
fitted to the spectrum that was taken as our template spec- 
trum in the cross-correlation (we have excluded the three 
Paschen lines that are possibly blended with absorption lines 
of the Ca II triplet). Hence, in order to convert the relative 
velocities in Fig.|^]and Table|2]to absolute velocities one has 
to add the systemic velocity. 

We have used the 32 I-band images that were obtained 
to acquire the source to construct a lightcurve. The images 
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Figure 3. Left panels: (a) The radial velocity curve derived by cross— correlating the Cir X— 1 spectra over the range 8700-8900 A (note 
that the velocities are relative to that of the template spectrum). Overplotted are the best— fit elliptical orbit treating the orbital period 
as a free parameter (light grey line) and that using the predicted orbital period from the X-ray dip ephemeris of Clarkson et al. (2004; 
dark grey line, P or f, =16.53 days; Note that our best— fit orbital period is consistent with that derived from the ephemeris of Clarkson 
et al. 2004). (b) The evolution of the /—band magnitude as derived from psf-fitting on the /-band acquisition images. We have plotted 
the relative errors in the magnitudes, in addition there is the ~ 0.2 magnitude uncertainty in the zeropoint (i.e. in the ordinate) (c) The 
contemporaneous X-ray flux as observed by the All Sky Monitor (ASM) on board the Rossi X-ray Timing Explorer. The points with 
horizontal error bars at an ASM count rate of 40 counts s _1 indicate the predicted time of the X— ray dip according to the X— ray dip 
ephemeris of Clarkson et al. (2004). Right panels: (d) The radial velocity curve folded on the X— ray dip orbital period of 16.54 days using 
the T from our best-fit to the radial velocity curve, (e) The phase folded I band magnitude, (f) The phase folded but unbinned ASM 
lightcurve using data 10 cycles before and after MJD 53475 (dwell data from MJD 53140—53807). The X— ray dip phase and T from the 
radial velocity curve are consistent at the 2 a level. 



were corrected for bias using the values from the overscan 
regions and flatfielded using sky flats taken within one or 
two days from the science images. To determine instrumen- 
tal magnitudes of Cir X-l and stars in its relative crowded 
vicinity, we used the point-s pread-function (psf) fitting rou- 
tines from DAOPHOT II dStetsorj|l987l) . running inside 
MIDAS. The instrumental magnitudes were placed in a com- 
mon photometric system by removing small magnitude dif- 
ferences between the different images, due to e.g. differences 
in the psf and the exposure times, by matching stars between 
these images and comparing their magnitudes. This common 
system was then calibrated using observations of standard 
star fields (PG 0942-029 an d L110; using the calibrated 
magnitudes bv IStetsonl 12000). which were imaged during 
photometric nights (March 27, April 11 and 12, 2005). As 
Cir X-l is only imaged in the /-band, we have not deter- 
mined colour terms. The uncertainty in the zeropoint is es- 
timated to be about 0.2 mag. The derived evolution of the 
/-band magnitude as a function of time and orbital phase 
is shown in the middle panels of Fig. It is clear that the 
/-band magnitude is constant except for a brief period that 
coincides with the X-ray dip and periastron passage. 



3 DISCUSSION 



We have obtained VLT phase resolved optical /-band spec- 
tra and images of the X-ray binary Cir X-l. The observed 
X-ray flux was low compared to that observed during other 



earlier spectroscopic observations (e.g. 


Johnston et alJll999l 


IJohnston et alJl200ll IClark et alJl200S 


) . If this lower X-ray 



flux corresponds to a lower intrinsic X-ray luminosity the 
optical light from the accretion disc thought to be caused 
in large by reprocessing of the X-ray luminosity might be 
lower as well. Hence, the putative companion star might 
be more readily observable now compared with high X-ray 
flux episodes. The spectrum we observed varies strongly as a 
function of orbital phase. Strong emission lines appear near 
periastron passage, filling-in the absorption line spectrum 
observable at other orbital phases. Cross-correlation of the 
absorption features in the 8700-8900 A range gives a ra- 
dial velocity curve with an orbital eccentricity of e = 0.45, 
a sin i — 16.9 lightseconds and a mass function of 0.019 
Mq. The question is; is this the radial velocity curve of 
the companion star, is it associated with the compact ob- 
ject/neutron star, or are the absorption lines caused by a 
circumbinary disc? 

The detected absorption line spectrum is consistent 
with a stellar spectrum if that star is a B5-A0 super- 
giant. This provides evide nce for a supergian t companion 
star in Cir X-l. Similarly. IClark et alJ i2003t) found emis- 
sion lines in the near-infrared spectrum of Cir X-l that 
often occur in mid-B supergia nts. Those results are con- 
sistent with our findings (but IClark et alJ 120031 favoured 
a scenario where the emiss ion features ar i se in the accre- 
tion disc/flow.) Previously, iMurdin et al] (^.980) proposed 
that the companion star of Cir X-l is an early type su- 
pergiant. However, the then known counterpart was later 
resolved into three stars red ucing the ma gnitude associ- 
ated with the companion star llMoneti lT992) which led peo- 
ple to discard the supergiant model for Cir X-l. On the 
other hand, an interstellar extinction of Ay = 5 was used. 
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Figure 2. From bottom to top: Spectra ranging from 8300— 
8900 A of Cir X— 1 (at ~phase 0.5; this work), the template star 
HD 33579 (spectral type A3Ia), HD 87737 (spectral type AOIb) 
and HD 77350 (spectral typ e A0III; the latter three spectra are 
from iLe Borgne et alj|2003) . If the absorption line spectrum is 
from the companion star of Cir X-l, the narrow Paschen lines 
show that the companion star must be a supergiant. The X-es 
mark regions where cosmic ray hits have been incompletely cor- 
rected for. 



Table 2. Relative radial velocities of the counterpart of Cir X— 1 
as displayed in Fig. [3^. To convert these velocities to absolute 
velocities add the systemic velocity of ~ —26 ± 3 km s — 1 . 



MJD 


Orbital phase 


Orbital phase 


Radial velocity 




(P orb = 16.53 d) 


(P orb = 16.68 d) 


(km s" 1 ) 11 


53445.3868 


0.323 


0.326 


-8.2±1.0 


53447.3749 


0.443 


0.446 


-18.1±1.0 


53448.2275 


0.495 


0.497 


-17.7±1.1 


53455.3966 


0.929 


0.926 


11.1±0.8 


53456.3680 


0.988 


0.985 


35.3±1.0 


53457.4003 


0.050 


0.047 


23.4±1.0 


53459.3877 


0.171 


0.166 


0.4±1.2 


53461.1724 


0.279 


0.273 


-8.9±1.1 


53465.3265 


0.530 


0.522 


-9.5±1.3 


53465.3517 


0.531 


0.523 


-6.9±1.4 


53465.3763 


0.533 


0.525 


-7.3±1.1 


53466.1605 


0.580 


0.572 


-12.4±1.1 


53471.2432 


0.888 


0.876 


22.0±0.8 


53472.2192 


0.947 


0.935 


30.4±0.8 


53473.2753 


0.011 


0.998 


49.7±1.1 


53491.1779 


0.094 


0.072 


27.2±1.0 


53492.2715 


0.160 


0.137 


8.1±1.1 


53493.2927 


0.222 


0.198 


11.5±1.6 


53495.2467 


0.340 


0.316 


-3.9±0.9 


53496.2661 


0.402 


0.377 


-4.3±1.0 


53498.3358 


0.527 


0.501 


-5.7±1.2 


53504.0955 


0.875 


0.846 


1.5±0.9 


53504.1153 


0.877 


0.847 


2.4±0.8 


53504.1351 


0.878 


0.848 


3.5±0.9 


53504.1549 


0.879 


0.850 


6.1±0.8 


53504.1748 


0.880 


0.851 


7.1±0.8 


53504.1946 


0.881 


0.852 


9.8±0.9 


53504.2144 


0.883 


0.853 


14.9±0.9 


53504.2343 


0.884 


0.854 


14.7±1.0 


53504.2541 


0.885 


0.856 


17.5±1.0 


53506.3040 


0.009 


0.978 


38.0±0.8 



a The error has to be increased by 6.5 to obtain a reduced \ 2 of 
Ril in the fit (see text). 



This value was derived on the basis of the lowest Nh mea- 
sured in the seventies. That Ng was determined from an 
X-ray spectrum obta ined from Aerobee r ocket data. No er- 
ror was given on Nh iMargon et al.lll97ll) . More recent N_s 
measurements from X-ray spectral fits with more sensitive 
satellites with a good soft response such as ROSAT, ASCA 
and Chandra give an N_g in the range 1.6-2 . 2X 10 22 cm" 2 
JPredehl fc Schmit,rlll99l iBrandt et alJll99ft ISchulzlll999t 
ISchulz fc BranddbOOa THowever. Ilaria et, alJ ^2003) again 
found a significant lower value for Nh modelling Beppo-SAX 
data. As we will show below there are three independant 
measurements that are consistent wi th the higher Nh fo r 
Cir X-l. First of all, as mentioned bv lBrandt et al.l |l996), 
a significantly lower value for Nh would not be compati- 
ble with the strength of the dust scatte ring halo in Cir X-l 
as found by IPredehl fc Schmittl {TllOV; . Secondly, the dis- 
tance ith^twouldbe derived from the low Nh value is ~4 
kpc lllaria et alJ 120051) . This distance is too low to explain 
the type I X-ray bursts, th ose require a lower limit to the 
distance of d=7.8-10.5 kpc Jjonker fc Nelemanslliool) . Fi- 
nally, the equivalent width of the diffuse interstellar band 
at 8620 A that we find in the optical spectrum is correlated 



Table 3. The best-fit orbital parameters. We provide both the 
solution for the fit with the orbital period as a free fit-parameter 
and that for where it was fixed to the orbital period extrapolated 
from the ephemeris of lClarkson et all 120041) . 



P (days) 
K (km s" 1 ) 

e 

w (deg) 

T (+53,475.0 days; MJD) 
a sin i (lightseconds) 
f(m) (M ) 



16.68±0.15 16.53 (fixed) 

25±2 
0.45±0.07 

2±12 
-1.7±0.4 
16.9±1.2 



25±2 
0.47±0.06 

-7±11 
-1.9±0.4 
16.8±2.0 
0.019±0.007 0.019±0.007 



with E(B_ V) l|Munarl2000h . This provides an A v = 12.2 for 
Cir X-l (taking R=3.1), consistent with the higher value de- 
rived from Nj? which usin g the conversion of Nh to Ay from 
IPredehl fc Schmitd dl995l) gives 9 < A v < 12. For these rea- 
sons and since the modellin g of the X-ray spe ctrum is model 
dependent as mentioned in llaria et alJ l)2005h we think that 
the higher value for Nh and hence Av is more likely to be 
right. 
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From the distance to Cir X-l from the type I X-ray 
bursts (d=7.8-10.5 kpc), Ay » 9-12, and 7=17.6, one gets 
an absolute magnitude —4.9 < My < —2.5 (taking into ac- 
count that V — I is nearly for mid-B /early A type stars, 
and using the relative extinction from ISchleeel et al]ll99Sl 
Ai — 0.6 xA^). An absolute magnitude of -2.5 would be too 
low and eve n -4.9 is on the low side for standard B5-A0 
supergiants fealona fc Cramptonll974T) . However, consider- 
ing the preceding binary evolution involving significant mass 
transfer from the neutron star progenitor to the initially less 
massive companion star and the current X-ray heating es- 
pecially at periastron the companion star of Cir X-l is not 
likely to be a standard supergiant. 

Let us investigate the consequences of the assumption 
that the absorption line spectrum and the radial velocity 
curve track the companion star of Cir X-l. Kepler's third 
law, the measured a op t sin i and P or 5 define a relation be- 
tween the inclination of the orbit with respect to the line- 
of-sight and the mass of the companion star (the solid lines 
in Fig. From the spectral classification B5-A0I we find 

2.l£ log g£ 2.4 JCoxl2000h. This gives = C \J^j^ 

(where C is in the range of 10-15). To avoid the neutron 
star going through the companion star at periastron, we 
constrain the radius of the companion star to be smaller 
than or equal to the periastron distance a (1 — e), with 
a sin i = —%r- — - a op t sin i and a op t sin i = 7.28 R . From 
Fig. Hthis yields M com p£l8.9M , £15M , £11.9 M 0) 
and <;1O.4M for M comp = 10, 5, 2 and 1.4 M , respec- 
tively. The inclination is constrained to be 7.1° <Jz<J 14.9°, 
8.9°£i£23.5°, 12.1°£i£50° and i£ 13.7° for a 10, 5, 2, 
and 1.4 M compact object, respectively. A normal B5-A0 
supergiant has a mass of ~10 M which would hence fit-in 
with a neutron star compact object of Cir X-l. On the other 
hand it is conceivable that if the companion star of Cir X-l 
has indeed spectral type B5-A0I it is no ordinary star due 
to the preceding evolution (as mentioned above). 

T he observed type I X-ray bursts (iTennant et alJ 
Il986bh are evidence of a neutron star nature of the com- 
pact object. Furthermore, as mentioned in the Introduction 
the neutron star has to have a magnetic field £ 10 12 Gauss. 
Also the presence of a strong radio jet is not compatible with 
a high magnetic field neutron star. These findings seem to 
be at odds with a supergiant (thus young) companion star. 
We conclude that if the presence of a supergiant companion 
star in Cir X-l is confirmed the neutron star magnetic field 
either decayed quickly to below £ 10 12 Gauss or the neutron 
star was born with such a low field. 

Assuming that a neutron star kick at birth did not give 
the binary system a large systemic velocity, one can calcu- 
late the radial velocity for Local Sta ndards of Rest along the 
direction of Cir X-l (see for instance lBerger fc Giesi200lfl to 
find that the systemic radial velocity of ~ —26 ±3 km s _1 at 
the location of Cir X-l gives a distance of either ~ 1.6 kpc or 
~11.8 kpc. A distance of 11.8 kpc is close to t he distance de- 
rived f rom the obse rved type I X-ray bursts llTennant et alJ 
Il986al: 7.8-10.5 kpc lJonker fc Nelemansll2004b . 

The optical I-band lightcurve (see Fig.[U) shows a clear 
brightening near periast ron. Th i s is si mila r in shape to the 
lightcurves published bv lMonetil (I1992T) and lGlassI (I1994T) . al- 
though one has to bear in mind that those lightcurves were 
obtained when the observed X-ray luminosity of the source 
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Figure 4. The mass of the companion star as a function of the 
orbital inclination with respect to the line— of-sight. The drawn 
lines are derived assuming that the compact object is from left to 
right, a 10, 5, 2, or a 1.4 M compact object. The dots are solu- 
tions if imposing the radius of the companion star to be smaller 
than or equal to the periastron distance a (1 — e) (we took the least 
constraining case of B5I, the radius of a A0I is somewhat larger 
hence the dots would all move to lower inclinations). For each of 
the solid line — dot combination the allowed parameter space is 
that on the solid line and below the dot. This yields the follow- 
ing constraints on the companion star mass for a 10, 5, 2 and 
1.4 M compact object; M comp ^ 18.9 M , £ 15 M , <;11.9M , 
and <J 10.4 M , respectively. The inclination is constrained to be 
7.1°^i^ 14.9°, 8.9°^i^23.5°, 12.1°^i<;50° and ip> 13.7° for a 
10, 5, 2, and 1.4 M compact object, respectively. Taking an AOI 
star instead of a B5I would change the inclination and compan- 
ion star mass constraint to 7.1°<Ji^ 10°, 7.3 M for the 10 M 
compact object. Under the assumption that the measured a sin i 
is not a op t sin i but ajvs sin *j the constraint on the companion 
star mass as a function of inclination is given by the dotted line 
(taking M NS = 1.4 M ). 



was much higher and that the data was phase folded using 
different ephemerides. The increase in I-band light corre- 
sponds to the phase where Paschen emission lines start to 
become apparent in the spectrum. Hence, it is probable that 
the enhanced mass transfer rate near periastron is respon- 
sible for both effects. For instance, the emission lines could 
be formed in the accretion stream from the companion star 
to the compact object and the enhanced I-band emission 
could be both due to these emission lines as well as due to 
enhanced continuum emission related to the stream impact 
site. 

In principle it cannot be ruled out that the absorp- 
tion line spectrum originates in the accretion disc. E.g. the 
/-band spectrum of the high inclination (accretion disc 
corona) source 2S 0921-630 shows Paschen absorption lines 
at orbital phases near 0.9 that could well be caused by 
the line of sight p assing through the accretion disc rim 
fl.Ionker et all2005r) . Those Paschen absorption lines are not 
present in the K1III spectral type that has been derived for 
the companion star in 2S 0921-630 (although X-ray heating 
effects may result in different spectral types being observed 
at different orbital phases, cf. t he observed spectr al type 
in Cyg X-2 changes from A5-F2 lCowlev et aJll97flh . How- 
ever, in radio observations of Cir X-l superluminal motion 
has been observed which limits the inclination of the jet axis 
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to the line of sight to i < 5° jFender et alJl2004l : the limit 
i < 5° has been derived assuming a distance of 6.5 kpc, if 
Cir X-l is indeed further away as is indicated by the burst 
properties, then the limit on i is more stringent still) . If our 
line-of-sight also goes through the accretion disc rim (at all 
orbital phases) it implies that the jet-axis is nearly in the 
plane of the orbit. Hence, it implies that the jet ploughs 
through the accretion disc if it originates close to the com- 
pact object. 

Nevertheless, if we assume that the a sin i we mea- 
sured is associated with the binary motion of the accretion 
disc/neutron star we can derive a limit on the mass of the 
companion star assuming the neutron star mass is 1.4 Mq. 
This limit is a function of the binary inclination and is given 
by the dotted line in Fig. 2] From the discussion above it 
would seem that the companion star mass is ~0.4 Mq since 
the binary inclination must be high in this scenario. A star 
of such a mass cannot have evolved off the main-sequence 
in a Hubble time unless the star was more massive initially 
(at least 0.8 Mq) and more than 0.4 Mq has been transfered 
already. However, it is unclear whether such an amount of 
mass can have been transfered while the orbital eccentricity 
is still 0.45. 

Finally, in this paragraph we consider other possible for- 
mation scenarios for the absorption line spectrum. E.g. in 
the peculiar X-ray binary S S 433 evidence for a circumbi- 
nary disc has been found jBhmdell et alJ 1200 J) , perhaps 
the absorption l i nes ar e formed in such a disc. However, 
iFilippenko et all (Il988t) have obtained I-band spectra of 
SS 433 and their spectra are markedly different from those 
that we observe from Cir X-l. In SS 433 the Paschen lines 
are double peaked and in emission whereas we find them 
to be in absorption except near periastron passage. On the 
other hand the inclinations at which we observe SS 433 and 
Cir X-l are differ ent. In this respect it is interesting to note 
that according to iTudose et alJ i2006T) we will be observing 
the system through the material that is swept up by the ap- 
proaching radio jet. However, it is difficult to imagine that 
that material is dense enough and has a velocity low enough 
to cause the narrow absorption lines. Finally, in all these 
scenario's the velocity changes of the absorption lines with 
orbital phase are difficult to explain. 

A potential way to determine whether the measured 
radial velocity curve is associated with the companion star 
is to obtain a high resolution spectrum. Such a spectrum 
could reveal narrow spectral features and allow for a better 
spectral classification e.g. the observed absorption feature 
near 8685 A consists of a blend of seve ral narrow lines if the 
spect rum is from the companion star jMunari fc Tomasellal 
1999) and if the companion star is an A-star the Ca II ab- 
sorption triplet can be separated from the Paschen lines. 
Furthermore, the rotational broadening of the spectral lines 
could be measured. Finally, due to precession of the orbit 
(i.e. apsidal motion), the anomalous period as measured by 
the time between periastron passages is different from the 
period that would be measured by eclipse timing (i.e. the 
sidereal period). Hence, if the X-ray dips in Cir X-l are 
due to (grazing) eclipses the sidereal period will over time 
differ from the anomalous period. 
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